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Quantum Chaos Versus Classical Chaos: Why is Quantum Chaos Weaker? 
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We discuss the questions: How to compare quantitatively classical chaos with quantum chaos? 
Which one is stronger? What are the underlying physical reasons? 
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; I I Chaotic systems occur in nature everywhere. Examples are: the turbulent flow of gases and liquids, e.g. 

I when pouring milk into a cup of coffee, in the bath tub, in hurricanes, or in the formation of galaxies. Chaos 
occurs in biology in the dynamics of populations and species. It occurs in the rate of heart beats, in attacks 
of epilepsy, and in the brain. It occurs also in chemical reactions. Most of the dynamics occuring in nature 
is chaotic. The occurence of non-chaotic systems in nature is more of an exception. The mathematical 
description of chaos has reached a mature state over the last decade. The tools being used are classical 
phase space, Lyapunov exponents, Poincare sections, Kolmogorov-Sinai entropy and others. 

There is also chaos in the quantum systems Q, |2; 01- People have investigated in quantum physics 
QQ . the analogues of classically chaotic systems. For example, a stadium-shaped billiard is classically chaotic. 
' The corresponding quantum system is a system of ultra-cold atoms bouncing against walls of stadium 
, shape, being created by interaction of the atoms with laser beams. Quantum chaos has been found to play 
' an important role in dynamical tunneling ^ |^ . The fingerprints of chaos in quantum systems and the 
mathematical tools of its description are quite different from those used in classical chaos. The neccessity 
for a different treatment is due to the nature of quantum mechanics: There is no proper phase space in 
' quantum systems - Heisenberg's uncertainty principle forbids that a point in phase space (uncertainty zero 
in position and momentum) exists. Heisenberg's uncertainty relation is a direct consequence of quantum 
mechanical fluctuations. A common approach to describe quantum chaos is random matrix theory and the 
use of energy level spacing distributions j^]. A fundamental conjecture by Bohigas, Giannoni and Schmit 
postulates that the energy level spacing distributions possesses a dominant part - which depends on the 
particular system - and a subleading universal part - independent of the particular system. The universal 
part gives a level spacing distribution of Wignerian type, if the corresponding classical system is fully 
chaotic. Such level spacing distribution can be generated also by random matrices of a certain symmetry 
class. The conjecture has not been rigorously proven yet, but has been verified and found to be valid in 
$-H ' almost all cases. 

This approach is successful and very popular. However, it has some short comings: First, the conjecture 
by Bohigas et al. holds strictly only in the case of a fully chaotic system, while in nature most systems 
are only partly chaotic, i.e. so-called mixed systems where its classical counter part has coexistence of 
regular and chaotic phase space. In such cases the quantum system yields a level spacing distribution, 
which is neither Wignerian nor Poissonian (where the latter corresponds to a completely regular system). 
There is no mathematical prediction of the functional form of such distribution. However, a number of 
interpolations between the Poisson and Wigner distribution have been proposed 0, By,aOlUJ- Second, 
we may ask: How about the comparison of the classical with the quantum system? And what is the 
quantitative degree of chaos? How can we answer this when the instruments used to measure chaos are 
quite different for both systems? 
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FIG. 1: Volume of chaotic phase space over total phase space versus energy. Comparison of results from classical 
action (dotted line) with quantum action (full line). 



Starting from this perspective and having in mind the goal to compare classical with quantum chaos, one 
may try the following strategy: Find a uniform description of chaotic phenomena, valid for both, classical 
and quantum systems. In more detail: Starting from non- linear dynamics and phase space in classical 
systems, one may look for a suitable quantum analogue phase space. Starting from random matrix theory 
and energy level spacing distributions of quantum systems, one may seek a random matrix description 
and a level spacing distribution of suitable dynamical objects in classical physics. In the following we will 
discuss some progress recently made in this direction. Using some of those results, we will compare for a 
particular system the chaotic behavior of the quantum system with the classsical system. The numerical 
analysis shows that the quantum system is globally less chaotic than the classical system. We believe that 
such finding is not limited to the particular system. In particular, we want to understand the underlying 
reason for such behavior. 



II. CASES WHERE QUANTUM CHAOS WAS FOUND TO BE WEAKER 

First, Casetti et al. considered the iV-component theory in the presence of an external field 
and in the limit of large N. They used mean field theory and observed a strong suppression of chaos in 
the quantum system, due to quantum corrections causing the system the system to move away from a 
hyperbolic fixed point rcsponsablc for classical chaos. Second, Matinyan and Miillcr studied another 
field theoretic model which is classically chaotic, namely massless scalar electrodynamics. They investigated 
the corresponding quantum field theory using effective field theory and loop expansion. They noticed that 
quantum corrections increase the threshold for chaos due to a modification of the ground state of the system. 
A third example is the kicked rotor which is a classically chaotic system in 1-D. Schwengelbeck and Faisal 
|14| considered the corresponding quantum system using Bohm's interpretation of quantum mechanics 
to introduce trajectories and a quantum equivalent phase space. They found that the Kolmogorov-Sinai 
entropy goes to zero in the quantum system, i.e. it is non-chaotic. This approach has been applied also to 
study chaos in anisotropic harmonic oscillators i coupled anharmonic oscillators and the hydrogen 
atom in an external electromagnetic field [13. Finally, Caron et al. [H considered the anharmonic 
oscillator in 2-D, which is classically a mixed chaotic system. Using the concept of the quantum action 
functional, a quantum analogue phase space has been constructed. As a result, the phase space portrait 
of the quantum system was found to be slightly but globally less chaotic for all energies (see Fig.pP). 



3 




-3-2-10 1 2 3 

X 



FIG. 2: Comparison of classical potential and quantum potential (multiplied by a scale factor 1.7331) for quartic 
coupling A = Also shown is the ground state wave function and two lowest energy levels. 
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FIG. 3: Comparison of classical instanton with instanton from the quantum action. 



Because the quantum action has been constructed from the classical action by taking into account quantum 
fluctuations [l^, 123, l2ll l2l, 12^ l2J, 12^ ■, hence the softening of chaos in the quantum system must be due 
to quantum fluctuations. We suspect that such softening effect may not solely show up in chaos. Indeed, 

looking at a double- well potential V(x) = -^(^^'^^ gj^ in the context of tunneling, we observe that quantum 

effects cause the quantum potential to be much "weaker" than the classical potential (see Fig.j5]), i.e. the 
potential wells are less pronounced and the potential barrier is much lower for the quantum potential (note 
that the quantum potential has a triple- well shape). The shape of the potential for tunneling translates 
into the shape of the instantons. Thus it comes as no surprise that the instanton of the quantum action 
(actually a double- instanton) is softer than the classical instanton (see Fig. [3]). 
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III. UNIFORM DESCRIPTION OF CHAOS 



Let us recall that the conjecture by Bohigas at al. is about random matrix theory and energy level spacing 
distributions to describe chaos in quantum systems, while chaos in classical systems is conventionally 
described in terms of phase space. That is, the tools to describe chaos in classical and quantum physics 
are different. For the purpose to better understand the physical content of this conjecture or eventually to 
find a mathematical proof, it is highly desirable to use the same language, respectively, tools in quantum 
physics as in classical physics. This can be viewed in two ways: Either one adopts the point of view 
that chaos should be analyzed in terms of phase space. Then a uniform description is achieved by use of 
classical phase space in classical physics and the quantum analogue phase space in quantum systems. How 
about the point of view that chaos should be analyzed in terms of random matrix theory and level spacing 
distributions? In the following we propose how to construct such a level spacing distribution, which will 
play the same role in classical physics as the energy level spacing distribution in the quantum system. 

One may immediately object that the energy level spacing distribution in quantum physics is due to the 
fact that one considers a system of bound states (with a " confining" potential) , which due to the rules of 
quantum mechanics gives a discrete spectrum. The discreteness of the spectrum is a quantum effect, i.e., a 
physical effect. Having in mind to construct a counterpart in classical physics, one may object that classical 
physics is continuous. There is no discreteness inherent in classical physics. If any discreteness occurs it 
will be due to some (mathematical) approximation. Hence, how can a level spacing distribution derived 
from a classical function be discrete and physical? First we want to propose a function and show how to 
construct a meaningful level spacing distribution. Afterwards, we will try to answer the last question. 

We would like to emphasize two properties: First, the level spacing distribution of the quantum system - 
in the case of a classically fully chaotic system - corresponds to a Wigner distribution. The type of Wigner 
distribution is determined by the symmetry of the Hamiltonian. E.g., there is the Gaussian orthogonal 
ensemble (GOE), the unitary ensemble and the symplectic ensemble. Thus, we need a function which 
carries the same symmetries. A function convenient for this purpose is the classical action. Another 
property to be emphasized is locality. Let us recall that the level spacing distribution, e.g. the GOE 
distribution, is invariant under orthogonal transformations. Because an orthogonal transformation maps 
any orthogonal basis of states onto another orthogonal basis, the level spacing distribution is essentially 
independent of the particular choice of the basis. That means, we may choose a basis that is quasi-local, 
i.e., built from square integrable functions which are identically zero everywhere except in a small interval 
where they are non zero and constant. Those box functions are almost local. The advantage of locality is 
the fact that symmetries of the Hamiltonian often have to do with transformations in position space, and 
locality facilitates the analysis of such symmetries. Why do we need to bother about symmetries? Because 
the energy level spacing distribution is meaningful only (and gives the Wigner distribution for a chaotic 
system, respectively a Poissonian distribution for an integrable system) if the energy levels have all the 
same quantum numbers (exept for the quantum number of energy). For example, for the spectrum of the 
hydrogen atom, one should take the bound states all with the same angular momentum quantum numbers, 
e.g., I ~ m = 0. 



We consider a system defined by a Langrange function L{q, q, t). Let S denote the corresponding action. 



Let qtraj denote the trajectory, i.e. the solution of the Euler-Lagrange equation of motion. Such trajectory 
is a function, which makes the action functional stationary. Each trajectory is specified by indicating the 
initial and final boundary points, i.e., q{t = 0) = gin and q{t = T) = qji. We assume that those boundary 
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points are taken from some finite set of nodes, Qin, Qfi £ {qi, ■ ■ ■ , Qn}- Thus with each pair of boundary 
points, {qk, qi) we associate a trajectory q^i'^'' it). Then we introduce an action matrix E, where the matrix 
element E^j corresponds to the value of the action S evaluated along the trajectory q%°'''{t)j 

E,,=5[C(t)]. (2) 

All matrix elements Efc; are real. They are also symmetric, E^; = Ejfc. Thus the action matrix E is a 
Hermitian N x N matrix. Consequently, E has a discrete spectrum of action eigenvalues, which are all 
real, 

cr(E) = {cti, . . .,aN} ■ (3) 



B. Symmetry 

For the purpose to compute a level spacing distribution from the action eigenvalues, one must first 
address the issue of symmetry. Let us consider for example the harmonic oscillator in 2-D. We may take 
the coordinates q^ to be located on the nodes of a regular grid (reaching from —A to +A on the x and y 
axis), with a spacing Aa; ~ Ay = a = const. For example, the action of the harmonic oscillator, with the 
classical trajectory going from Xa to Xb in time T is given by 

^a,b= J^^^[ixl+xl)cOsiLuT)-2Xa-Xk] . (4) 

This function is invariant under rotations. When choosing the coordinates qk to be located on the regular 
grid, the continuous symmetry of rotation will become a discrete symmetry of finite rotations (a group). 
We have to find the irreducible representations of such group and sort the action eigenvalues according 
to those irreducible representations (this can be done by inspecting the properties of the corresponding 
eigenvector). Then one has to select a particular representation and retain a subset of eigenvalues in that 
representation. The action level spacing distribution can then be obtained from the action eigenvalues in 
such subset. 

This procedure is feasible. However, it has two disadvantages. First, finding the irreducible represen- 
tations and classifying the eigenvectors accordingly is laborious. More importantly, the fact that one has 
to work with a subset of eigcnvualcs only and not the whole ensemble of eigenvalues means a drastic re- 
duction of the size of the statistical ensemble. In other words, the statistics of the resulting level spacing 
distribution will deteriorate. 

For those reasons it would be highly desirable to avoid the above strategy. This is indeed possible by 
using the following trick. One can camouflage the symmetry by choosing the coordinates off the nodes of 
the regular grid. That means, for example to define new coordinates as follows, 

deform , /r-\ 

qk =qk + efc , (5) 

where e^. denotes a randomly chosen small deformation (in angle and length). In this way the nodes are 
irregularly distributed. Consequently, the discrete symmetry of the action matrix E^j disappears when 
replacing it by the "deformed" action matrix, 

j^'l^'^f"™ ~ action evaluated along the classical trajectory 

r J deform , j deform //^\ 

irom node q^, to node q^ . (o) 



In doing so we avoid a laborious symmetry analysis and secondly will have a better statistics! 
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C. Action level spacing distribution 



Now wc want to construct a level spacing distribution of action levels. We proceed in analogy to random 
matrix theory and the method of constructing an energy level spacing distribution. For an overview on how 
to compute level spacing distributions see Ref.[2^. One has to separate the dominant system dependent 
part from the subleading universal part which describes the properly normalized fluctuations. Because 
we are here interested only in the fluctuation part, we suppress the leading part. One should note that 
this means to discard all physical information which depends on the particular system. For example 
thermodynamical functions can not be computed from the subleading fluctuating part. The strategy to 
obtain the subleading part is called unfolding. One constructs a fit to the original spectrum and multiplies 
the spectrum such that on average the mean spacing distribution becomes unity. Also the integrated level 
spacing distribution will be normalized to unity. 

We have applied this to simple integrable systems in 1-D. For integrable systems one would expext a 
Poissonian distribution for the action level spacing distribution. Preliminary results are compatible with 
a Poissonian distribution. The following remarks are in order. The first numerical results for integrable 
systems have to be repeated with precision and analyzed carefully. Second, one wants to see what happens 
in chaotic systems. Possibly such strategy applied to a fully chaotic system will result in a Wignerian action 
level spacing distribution. Finding an answer will be computationally much more involved, simply because 
the action functions for the integrable systems considered above are analytically known, while for a chaotic 
system (i.e. non-integrable) this needs to be calculated numerically. Moreover, the numerical precision 
required needs to be sufficient to resolve small fluctuations. In the statistical sense, one is interested in 
a sample of large size. But that means that after unfolding the fluctuations will become small and hence 
require a high numerical precision for its resolution. Presently, numerical studies of such question are under 
way. 

Let us get back to the question posed at the beginning of Sect lllll How can a level spacing distribution 
derived from a classical function be discrete and physical? In our opinion, the answer lies in the fact that 
the level spacing distribution is universal, that is, it does not depend, for example, on the parameters of 
the discrete grid. Different grids give the same result. This has been verified numerically. It also should 
not depend on the deformation (as long it is not too close to zero and the discrete symmetry is restored). 
Also this has been verified numerically and found to be satisfied. Thus one could in principle go with the 
volume of the lattice V = (2A)^ to infinity and with the lattice spacing Ax = Ay to zero. The result should 
not change, but one would have reached the continuum limit. The discreteness would then disappear. 

IV. RENORMALISATION FLOW OF PARAMETERS OF THE QUANTUM ACTION 

In Sect^we have seen examples for the observation that quantum chaos seems to be weaker than chaos 
in the corresponding classical system. Of course it would be interesting to explore a much wider class of 
systems in order to see if such observation holds more generally. Here we want to pick one of the above 
examples, namely the the chaotic anharmonic oscillator in 2-D and try to understand why quantum chaos 
is weaker than classical chaos. The classical action is given 



For A = the system is reduced to the standard harmonic oscillator, which is integrable. The chaoticity 
is introduced and controlled by the parameter A. Thus small A causes mild chaos, while large A makes the 
system strongly chaotic. The quantum action has been postulated to be of the functional form like the 
classical action, i.e. the kinetic term of the quantum action may differ in the value of the mass, and the 





,2 



V 



(7) 



7 



potential of the quantum action should also be local, depend only on coordinates, but may have a different 
functional form. Here let us consider an ansatz of the following form 

S = dt^[x^ + f)-V{x,y) , 

V = Vq + V2{x^ + y^) + V22x'^y^ + higher order polynomials . (8) 

As a quantitative measure of the strength of chaos we take the strength of the parameters of the action, in 
particular, the parameter A. The parameters of the quantum action can be interpreted as a " renormalisation 
effect" of the parameters of the classical action. The calculation of those parameters has to be done 
numerically, following the definition of the quantum action to be a functional which fits the transition 
amplitudes |l9l| . However, in a certain limit, the quantum action is known to be an exact representation of 
the transition amplitudes and moreover the action is related via differential equations to the classical action 
p2ij . This limiting case is using imaginary time and let time go to infinity (Feynman-Kac limit). Because 
we want to obtain an analytical result, we will use perturbation theory. This means that we consider the 
regime of small A = V22- 

In order to simplify the matter, let us start by considering the system in 1-D. Thus we have the potential 

V{x) — —muj'^x'^ + Ax'* = V2x'^ + v^x'^ . (9) 

We assume A to be small, 

— - « 1 , 10 

V2 

where A^c introduces a physical length scale, e.g. the analogue of the Bohr radius. According to Ref . [2^ 
the following relation between the classical and the quantum potential holds, 

2miV{x) - Egr) = 2MVix) - io) - ^i^^^M=M sgn(x) . (11) 

^2m(T/(x) - vo) 

We define the functions 

W{x) ^ 2m{V{x) - Egr) , (12) 
U{x) = 2m{V{x) - Wo) , (13) 

where Egr denotes the ground state energy of the lowest eigen state of the quantum mechanical system 
(obtained from the Schrodinger equation). Due to Eq.lQ, the function W{x) must have the form 

W{x) = Wo + W2X^ + wj^x'^ . (14) 

This function is symmetric with respect to parity. Then the function U {x) representing the quantum 
potential, will be parity symmetric also. We make an Ansatz of the form 

U (x) = Uq + U2X^ + Uix"^ + UqX^ + • • • (15) 

The assumption that the expansion parameter A is small is now expressed by 

W4 = wf^e, £<<! (16) 



Now using Eas. H14ll5|) in combination with Ea. lll|l . we obtain 

+ UqX^ 



2 + Uix"^ + uga;'* + • • • 



for .T > . (17) 



The smallness of W4 imphes that the terms of fourth order and higher in x occuring in the function U {x) 
are small compared to the second order, i.e. 

UiX^ + Uqx'^ + • • • «U2 ■ (18) 

Now doing a Taylor expansion in the small terms u^^x"^ + ugx* + ■ ■ • allows to express the r.h.s. of Ea. p7|l 
as a polynomial in x. Then comparing terms in x order by order, we find the following relations 

^ A 2 
4^/?I2 



h 

W2 = U2 T 6U2U4 , 

4^M 



W4 = U4 3 (IOM2M6 — 4^4) . (19) 



Now we try to find the parameters U2 and M4 as solution of those equations. The first equation gives 

U2 = (wo/hf . (20) 

On the other hand we have, due to Ea. H13|) . 

Wo — —2mEgr ■ (21) 

Due to the anharmonic perturbation the ground state energy is different from the ground state energy 
of the harmonic oscillator. However, because of the smallness of the perturbation, we can express the 
groundstate energy Egr using perturbation theory as a power series in e, 

Egr = + eE^^^ + £2^(2) ^ . . . ^ (22) 
where E'-^'^ = E^-^". Thus from Eas. (|21l22|) we obtain 

2mE^^^ 

^muj + e + O(e^) , (23) 

or 

U2 = m^bj^ + e V 0(e 

h 



2^(1 

W2[l + 

Next let us consider Ea. p9b '). We obtain 

' ~ " 3ft 



^«2[l + -pJ^e + 0(e2)] . (24) 



2 

U4 = --zT\fu^{w2 - U2) ■ (25) 



Recalling W2 = m^w^ and Ea. lf^ . we find 

2 777,3^3 £;(!) 



In Eas. H24l26|l we have expressed the parameters of the quantum action in terms of the parameters of the 
classical action. However, it remains to compute the energy E^^K Again we use stationary perturbation 
theory. The Hamiltonian is given by, taking into account Eas. (|14ll6|l . 



H 
Hio) 

Hii) 



H 



(0) 



(1) 



F , ^ 2 2 

— muj X 

2m 2 

(0) 



2m 



To first order of perturbation theory in e the energy E^^'> is given by 



which yields the result 



(0) 



Substituting this result into Eas. (|24l26|l wc finally obtain 

3^1 



U2 = W2 



1 + 



)wi + Oie') 



and 



1 



U4 



(27) 



(28) 



(29) 



(30) 



(31) 



V. INTERPRETATION 

Let us see what happens when we keep the classical parameters fixed, exept for u'4, i.e., we keep w'^^'^ 

fixed and vary e. Note that ^2 > and W4, > 0. We also have w'^^ > and e > 0. Now we want to study 
what happens when e ^ 0. Ea. (|30|l yields 

U2 > W2 
U2 >e~>0 W2 

Likewise, Ea. l(^ yields 

Ui >e->0 Wi ^£-^0 . (33) 

In other words, in the limit e ^ the classical potential approaches the potential of the harmonic oscillator. 
The potential of the quantum action asymptotically approaches the classical potential, hence also the 
harmonic oscillator potential. That is, the renormalisation group fiow of the parameters U2(e),M4(e) goes 
to a Gaussian fixed point. Second, for any value of e the value the quadratic term of the potential is larger 
for the quantum potential than for the classical potential. Third, for any value of e the value the quartic 
term of the potential is smaller for the quantum potential than for the classical potential. Recall that the 
quadratic term is the term, which, if it would stand alone, would make the system integrable. On the other 



(32) 
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hand, the quartic term is the term which drives the system away from integrabihty (and introduces chaos 
in 2-D). Thus we find that quantum fluctuations, which are the cause for the differences A2 = W2 — U2 
and A4 = W4 — U4 to be non-zero, have the tendency to drive the quantum system closer to the regime of 
integrabihty. 

Now the above perturbative calculations were performed in 1-D. Chaos in time-independent Hamilton 
systems exists only for D > 2. A similar, but more tedious calculation can be performed in 2-D. ft confirms 
the above result that quantum fiuctuations drive the system closer to the regime of integrabihty and away 
from the regime of chaos. 

Of course such perturbative calculations are meaningful only in a neighborhood of the Gaussian fixed 
point. It would be desirable to extend the calculations to a larger regime. However going to higher order 
of perturbation theory would make those calculation much more tedious. Nevertheless the perturbative 
result gives some insight into the dynamical consequences of quantum fluctuations. 
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